Effects of tensor interaction on pseudospin energy splitting and shell correction 
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In the framework of a Skyrme-Hartree-Fock approach combined with BCS method, the role of the 
tensor force on the pseudospin energy sphtting for tin isotope chain is investigated. The tensor force 
turns out to obviously affect the pseudospin energy splitting of the spin-unsaturated nuclei. Since 
the tensor force shifts the single-particle levels, it modifies the single-particle level density and the 
shell correction energy thereof. The influence of the tensor interaction on shell correction energy is 
considerable according to our analysis taking a magic nucleus *^^Sn as well as a superheavy nucleus 
■^^^114 as examples. This modification of the shell correction energy due to the tensor component 
affects the stability of the superheavy nuclei. 
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I. INTRODUCTION 

Properties of exotic nuclei are not only interesting from 
the viewpoint of nuclear structure, but also have impor- 
tant implications in nuclear astrophysics, such as in nova 
and supernova explosions. X-ray bursts associated with 
explosive hydrogen burning, rapid proton capture pro- 
cesses, etc. In view of the recent rapid progress in the 
discovery of exotic nuclei driven by the advent of the new 
generation of radioactive ion beam facilities, an impor- 
tant problem is to understand how the shell structure 
evolves from stable to exotic nuclei, in which the role of 
the tensor force has received wide attention. The effect 
of the tensor force on the shell evolution has been ex- 
plored by Otsuka and his collaborators within the shell 
model. They found that the tensor force strongly affects 
the evolution of the single-particle energy level spacing 
and shell structure 

In the case of Skyrme interactions, a zero-range ten- 
sor component was actually present in the original pa- 
pers of Skyrme However, this component was omit- 
ted when realistic Skyrme parameter sets were deter- 
mined, and it has been systematically neglected for sev- 
eral decades in practical applications. Until very recently, 
the tensor component was added to the existing Skyrme 
parametrizations or included by fitting new parametriza- 
tions to investigate the contributions of the tensor force 
to the spin-orbit splitting evolution of single- 

particle energies , shell evolution d, H, |^, stability 
of superheavy nucl ei [ill, hinging energy [l2l . [l3j , nu- 
clear deformability [l3l Il4l|. properties of excited states 
[Tsl [TgI , and the spin and spin-isospin instabilities of nu- 
clear matter [l7| . 
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The concept of pseudospin was supported by experi- 
mental observations about 40 years ago that the single- 
nucleon doublet levels with quantum numbers {rir, I, 
j = 1 + 1/2) and (n^ — 1, I + 2, j = I + 3/2) lie very close 
in energy, where n^, I, j are the single nucleon radial, 
orbital, and total angular momentum quantum numbers, 
respectively. This single-nucleon doublet pair can be re- 
labelled as a pseudospin doublet: (n^ — rir — l, I — l + l, 
j — I ±1/2). Then the two states in the doublet are al- 
most degenerate with respect to pseudospin. This sym- 
metry has been used to explain a number of phenomena 
in nuclear structure including the deformation [l8l [. su- 
perdeformation [3, magnetic moment [13, [lH and iden- 
tical bands |22h24|. It is also pointed out that the con- 
servation of the pseudospin symmetry plays an essential 
role in stabilizing the neutron halo structure and gener- 
ating the neutron shell effects when the neutron num- 
ber approaches the neutron drip line [25|. The details 
about the pseudospin symmetry can be found in Refs. 
[26l - l28j . It is interesting to discuss the influence of the 
tensor force on the pseudospin symmetry. In this work, 
we shall discuss the effect of the tensor component of the 
two-body effective interaction on the pseudospin energy 
splitting by employing a Skyrme-Hartree-Fock approach, 
taking a long chain of Sn isotopes as an example. In ad- 
dition, the effect of tensor interaction on shell correction 
energy will be investigated. The liquid drop model ex- 
plains well the gross features of nuclear fission and mass, 
however, there are systematic deviations from the smooth 
mass. The shell correction, which is a fluctuation in the 
binding energy, can be supplemented to the liquid-drop 
model to improve the description of the nuclear masses 
and fission barriers 29-35]. For superheavy nuclei, as an 
important property, the shell correction energy is widely 
investigated because it is relevant to the stability of a su- 
perheavy nucleus. The strutinsky method [s^, which is 
widely used today, will be employed to extract the shell 
correction energy here. With the single-particle spectra 
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obtained from the Skyrme-Hartree-Fock approach, the 
shell correction energies for heavy and superheavy nuclei 
can be calculated and one can analyze the effect of the 
tensor force. 

This paper is organized as follows: In section II, a brief 
theoretical method is presented. The effects of the tensor 
interaction on the pseudospin energy splitting and shell 
correction are discussed in sections III and IV, respec- 
tively. Finally, a brief summary is provided in section 
V. 



II. TENSOR FORCE IN THE FRAMEWORK OF 
THE SKYRME INTERACTION 



The tensorpart of the Skyrme effective interaction is 
written as ^ 



VT = ^ \{<Ji-k'){a2-k') 



(cTi • k)(o-2 • k) - i(CTi • cr2)k^ 



T 

-C/[(ai-k')(5(ri-r2)(a2-k)] 
-^(ai.a2)[k'.<5(ri-r2)k]. 
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1-2 
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(1) 



where k = ('^1 — '^2)/(2i) acts on the right and k' — 
2) I {"2.1) acts on the left. T and U provide 
the intensity of the tensor force in even and odd states 
of relative motion, respectively. The spin-orbit potential 
including the tensor contributions is given by jSj 



2r \ dr dr 



J-H+P^-^], (2) 
r r 



where q{q') denotes the like (unlike) particles and the Jq 
is the proton or neutron spin-orbit density given as 



(3) 

where is the BCS occupation probability of each orbit. 
The orbital with j> = Z-l-1/2 gives a positive contribution 
to Jq while the orbital with j< =1 — 1/2 gives a negative 
contribution to Jq. The strengths of these contributions 
a and /3 are expressed as a = ac -\- ctT and 13 = (3c + (3t- 
ac and /3c are related to the central exchange part of 
the interaction and can be written in terms of the usual 
Skyrme parameters 
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I3c = --^(iia^i + ^2a;2). 



(4) 



ax and /3t are related to the tensor part in the following 
way 
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(T + U). 



(5) 



The central exchange and tensor contributions modify 
the energy density H by 



AH =-a- 
2 



Jl + Jl)+pJnJp- 



(6) 



One set of the parameters employed in the present work 
is the same as that in Ref. Q: ay = —170 MeV fm^ 
and = 100 MeV fm^ added perturbatively to the 
Skyrme force SLy5 iSTj (marked by SLy5-|-TF), which 
can fairly well explain the isospin-dependence of the ex- 
perimental energy differences between the single-proton 
states outside the Z = 50 core for Sn isotopes, and the 
single-neutron states outside the Z = 82 core for Z = 82 
isotones both quantitatively and qualitatively. Recently, 
Zou et al. analyzed the evolution of the spin-orbit split- 
tings in the Ca isotopes and A*" = 28 isotones with this 
SLy5-|-TF and they found that adding the tensor contri- 
bution can qualitatively explain in most cases the empir- 
ical trends (7|]. 
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FIG. 1: Energy differences between the I/111/2 and lg7/2 
single-proton states along the Sn isotopes. The experimental 
data are taken from Ref. [S^. For the proton-rich isotopes, 
one of or both of these two levels lie on the positive energy 
states (unbound), so we do not present them here. 



III. EFFECT OF THE TENSOR FORCE ON 
THE PSEUDOSPIN ENERGY SPLITTING 



Apart from affecting the spin symmetry, the tensor 
force influences the pseudospin symmetry. The shift 
of single-particle levels might have different origins and 
mechanisms. In the present study, we focus on the ef- 
fect of the tensor force. The calculations are performed 
in a coordinate space using a box size of 20 fm and a 
mesh size 0.1 fm. The pairing correlation is accounted 
in the BCS formalism with an energy gap determined 
by the observed odd-even mass differences for an open 
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FIG. 2: Proton pseudospin orbit splitting Ae measured by 
e(2p3/2) -£(1/5/2) for Idand by £(137/2) -e(2d5/2) for 1/in 
the Sn isotopes. 



shell and being zero for a closed shell. In addition to the 
SLy5+TF, the parameter set TBI including the tensor 
force (T31+TF) is employed in the present study, which 
also explains the experimental energy differences between 
the I/111/2 and 1177/2 single-proton states along the Sn iso- 
topes, while the interaction T31 without the tensor com- 
ponent fails to reproduce the trend of the experimental 
data, as shown in Fig. 1. Note that the calculations with 
the SLy5-|-TF were firstly performed in Ref. Q using 
other pairing formalisms, yet their results coincide with 
our calculated ones, which is therefore consistent with 
the conclusion that the pairing effect is negligible for the 
spin-orbit splitting suggested in Ref. [s^ . We display the 
proton pseudospin energy splitting of Id (2^3/2,1/5/2) 

and 1/ (I57/2, 2^5/2) of the Sn isotopes in Fig. 2. The re- 
sults with the parameter sets SLy5 (SLy5-|-TF) and T31 
(T31-I-TF) show the same trend. The proton spin current 
does not contribute to isospin dependence of the spin- 
orbit potential for an isotope chain, however, the neu- 
tron spin current is uniquely responsible for the isospin 
dependence and affected greatly by the sign of the pa- 
rameter /3t. From TV = 56 to TV = 64, J„ is reduced as 
the 1(77/2 neutron orbit is gradually filled. Because of a 
positive value of /?t, the absolute values of the proton 
spin-orbit splittings are enlarged. Therefore, the energy 
level of the proton orbit I/5/2 is pushed up while that 
of the 2p3/2 is pulled down, leading to the reduction of 

the pseudospin orbit splitting of Id. When the neutron 
number iV changes from 66 to 70, the filling of the 2c?3/2 



FIG. 3: Same as Fig. 1, but for neutrons. 



neutron orbit also reduces the pseudospin orbit splitting 
of Id. The pseudospin energy splitting of 1/ should in- 
crease in this region, nevertheless, the lower panel of Fig. 
2 shows an opposite trend. The reason for this result 
is that the centroid of the \g levels shifts down more 
rapidly than that of the 2d levels as the neutron number 
increases, which submerges the weak isospin-dependence 
of the tensor effect in this region. Moreover, from N — 72 
to N = 82, the I/111/2 orbit is gradually filled, which re- 
duces the spin-orbit splitting, so that the I/5/2 and 1177/2 
orbits shift downwards, but the 2p3/2 and 2d5/2 shift up- 
wards. Accordingly, the pseudospin energy splitting of 
Id is enhanced, however, that of the If is reduced as the 
neutron number increases. 

Although the proton spin current does not affect 
isospin dependence of the spin-orbit potential for a given 
isotope chain, being spin-unsaturated for protons, the 
filling of the proton orbits affects nucleonic spin-orbit 
splitting as compared to that without the tensor force. 
The negative value of ax and positive value of Jp en- 
hance the absolute value of the spin-orbit potential and 
hence increase the spin-orbit splitting. From = 56 to 
A^ = 68, the tensor effects resulting from the proton spin 
current (ax term) and neutron spin current (/Sx term) 
affect the spin-orbit splitting in the same way. Conse- 
quently, the pseudospin energy splitting of Id is reduced 
and 1/ is enlarged evidently compared to those without 
the tensor force. From A^ = 70 to A^ = 82, the ar and 
terms vary in an opposite way, and the contribution 
of the Ut term is cancelled out gradually. Accordingly, 
the difference of the pseudospin energy splitting between 
that with and without the tensor force is reduced gradu- 
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ally. On the whole, compared to that without the tensor 
force, the pseudospin symmetry is recovered to a large 
extent for the pseudospin doublet of Id but is broken for 
If. The tensor interaction thus brings a distinct change 
in pseudospin symmetry. For neutrons, the role of the 
tensor force is more complicated due to the variation of 
neutron number, but the obvious alternations caused by 
the tensor component are distinguished, as shown in Fig. 
3. The tensor force leads to the energy inversion of some 
pseudospin doublets, such as (2p3/2, 1/5/2) when N varies 
from = 56 to 60 as can be seen in Fig. 3. 

For the pseudospin doublets except 11^3/2)1 

such as (1^7/2, I59/2) and (1/19/2, lft-11/2) in the spin- 
unsaturated heavy nucleus, they are all composed of a 
j> =1 + 1/2 orbit and a j < =1 — 1/2 orbit. Accord- 
ingly, their pseudospin symmetry could be affected con- 
siderably by the tensor effect because the two orbits are 
shifted by the tensor component in an opposite direc- 
tion. For the pseudospin doublets {npi/2, ^P3/2), only 
the np3/2 orbit is shifted due to the tensor force (the 
npi/2 orbit is not affected due to its zero orbit angu- 
lar momentum), and hence the tensor force also affects 
these pseudospin splittings. Our aim is to illustrate the 
significant influence of the tensor force on the pseudospin 
symmetry. Because the pseudospin symmetry is affected 
by the tensor force, a number of phenomena [T8l - [25j | be- 
ing related to this symmetry can be accordingly affected. 
In addition, the rotation bands of the deformed nuclei 
based on the single-particle levels should be influenced 
by the tensor force, which needs investigation but it is 
not going to be easy because of the breaking of the time 
reversal symmetry. 

IV. EFFECT OF THE TENSOR FORCE ON 
THE SHELL CORRECTION 

Once the single-particle spectra are obtained, the shell 
correction energies can be calculated. Calculating such 
a correction requires the subtraction of the "smoothing 
varying" part of the sum of single-particle energies 

N{Z) 1 

EsheW = E ~ E = ^ eiHi - / eg{e)de, (7) 
1=1 

where N{Z) is the neutron (proton) number, and are 
the single-particle energies. The neutrons and protons 
are treated quite separately. The smoothed Fermi level 

A is obtained from the equation N{Z) = J^^g{e)de. 
Taking into account the n^-fold degeneracies of the single 
particle orbits, the smoothed level density g{e) takes the 
form of 

g{e) = If^nj f ^) , (8) 

where 7 is the smoothing range in units of Hujq = 41^"^/'^ 
MeV. The folding function is usually taken as f{x) = 



-^e ^'^ {x^) , where l)I,^{x'^) is an associated La- 
guerre polynomial. Here the energy cutoff of the single- 
particle spectra is selected to be 40 MeV. 

It remains to determine the order of the associated La- 
guerre polynomial M and the smoothing range 7. Ideally, 
the calculated EsheW should be not dependent on the spe- 
cific values in a broad range of reasonable values, namely, 
the so-called plateau condition 

In our calculations, it is found that A/ = 3 is an optimal 
value just like that in Ref. [l^l, in which the associ- 
ated Laguerre polynomial takes the form of Ly^{x^) = 
-x^/Q + lx^/\ ~ 35x78 -f 35/16. Fig. 4 presents the 
neutron and proton shell correction energies as a func- 
tion of the smoothing range 7 taking the nucleus ^"^^Sn 
as an example, where the Skyrme-Hartree-Fock equations 
are solved with box sizes of i?box = 8 ~ 12 fni with 
the parameter sets SLy5 and T31. One can find that 
the plateau condition is well fulfilled within the range 
1.5 ^ 7 ^ 2.7 for neutrons and 1.2 ^ 7 ^ 2.0 for protons 
in the case of i?box = 10 fm. The plateau vanishes gradu- 
ally as i?box > 10 fm. And i?box = 10 fm is large enough 
in our calculations to achieve the convergence for the to- 
tal binding energy with an accuracy of 0.03%. When 
-Rbox — 10 fm, the calculated shell correction energies for 
^■^^Sn with the parameter sets SLy5 and T31 in the cases 
with and without the tensor force, are presented in the 
insets of Fig. 4. This shell correction shows a parameter 
dependence. Without the tensor force, the shell correc- 
tion energies for neutron (proton) are about -12.0 MeV 
(-7.0 MeV) and -10.3 MeV (-5.5 MeV) by adopting the 
parameter sets SLy5 and T31, respectively. The tensor 
force shifts the single-particle levels, thereby the single- 
particle level density undergoes a corresponding change, 
so that the shell correction energy is modified. It can be 
found that the shell correction energies for both neutrons 
and protons are reduced by the tensor force, the reason 
for which is partly because the tensor interaction makes 
the energy level density near the Fermi surface become 
lower. In order to show a clearer change in the shell cor- 
rection energy caused by the tensor force, we plot in Fig. 
5 the net effect of the tensor force on the shell correc- 
tion energy Ai?, namely, the difference between the shell 
correction energies with and without the tensor force. 
As can be seen, the Ai? presents a more broad range of 
7 from 1.5 to 3.5 to fulfill the "plateau condition" with 
-Rbox = 10 fm. With i?box = 11 fm, though the neu- 
tron and proton shell correction energies do not show the 
"plateau "as presented in Fig. 4, the net contribution of 
the tensor force to the shell correction energies shows the 
"plateau " , especially for protons and nearly provides the 
same result as that with 7?box = 10 fm. This indicates 
the reliability of our calculations to a large extent. The 
error caused by the box size can be cancelled out to a 
large extent when one computes the contributions from 
the tensor force. One can find that due to the tensor 
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FIG. 4: Neutron and proton shell correction energies as a function of the smoothing range 7 for Sn with different box sizes. 
The insets display the shell correction energies with and without the tensor component with the box size -Rbox = 10 fm and 
their coordinate axes are the same as the figures. 
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FIG. 5: Modification of neutron and proton shell correction 
energies due to the tensor force as a function of the smoothing 
range 7 for ^^^Sn. 



force, the neutron shell correction energy is lowered by 
about 2.3 MeV and 3.0 MeV, and the proton shell cor- 
rection energy is lowered by about 1.1 MeV and 1.8 MeV 
with the parameter sets SLy5+TF and T31-I-TF, respec- 
tively. In other words, the effect of the tensor force on 
the shell correction is considerable. By the way, the shell 
correction energy at the saddle point turns out to be too 
important to be neglected 41], so the tensor force may 
play a non-negligible role in fission processes. 

Now we turn to the shell correction for superheavy 
nuclei. The existence of the superheavy island was pre- 
dicted in 1960s. However, there is no consensus among 
theorists with regard to what should be the next doubly 
magic nucleus beyond ^°^Pb. Nearly all of modern calcu- 
lations predict the existence of a closed neutron shell at 
N = 184. However, they differ in predicting the position 
of the closed proton shell. Modern calculations predict 
that they appear at Z = 114, 120, 124, or 126 [s^, IH- 
li?! depending on the models employed. We investigate 
the tensor effect on the shell correction of superheavy 
nuclei taking 298114 as an example. The influence of the 
tensor force on the neutron and proton shell correction 
energies (i?box — 12.5 fm) is illustrated in Fig. 6. With 
the parameter sets used above, we find that the shell cor- 
rection energy for protons is altered by 4 ~ 5 MeV due 
to the inclusion of the tensor component, and is much 
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FIG. 6: Modification of neutron and proton shell correction 
energies due to the tensor force as a function of the smoothing 



range 7 for ^"llA. 



larger than the 1.0 ^ 1.6 MeV for neutrons. This im- 
phes that 2^^114 becomes more stable due to the pres- 
ence of the tensor force with the present parameter sets 
since the shell correction energy can reflect the stability 
of nuclei. We plot in Fig. 7 the single-particle levels 
near the Fermi surface to show the effect of the tensor 
force (i?box = 20 fm). As can be seen, the shell gap at 
Z = 114 is obviously enlarged due to the tensor force, 
and the shell gap at = 184 is changed slightly, which 
is consistent with the above discussion about the shell 
correction energy and is agreement with the results in 
Ref. IJJJ. The difference between the 4si/2 and 2/iii/2 
neutron orbits determines the size of the gap at N = 184, 
yet the tensor force only directly affects the 2/iii/2 orbit. 
In consideration of the fact that the ax and /3t terms in 
Eq. ([2]) take effect in an opposite way, their contributions 
to the spin-orbit splitting are cancelled out to a certain 
extent, hence the predicted neutron shell structures is 
modified only slightly by inclusion of the tensor compo- 
nent. However, for Z = 114, the two proton orbits 2/5/2 
and 1*13/2 shift in an opposite direction due to the ten- 
sor component, thus the gap is changed obviously. The 
discussions here aim at displaying the effect of the ten- 
sor force on the structure and stability of the superheavy 
nucleus. The theoretical results at least imply that the 
tensor force should play an important role in the shell 
structure of the superheavy nucleus. We would like to 
say that many-body approaches with various parameter 
sets tend to predict different locations of the shell gaps, 
which perhaps results from the uncompleted knowledge 
about the tensor force to a large extent. 



-2 F 



4 - 2f, 



-6 - 



> 



S -8 - Ih 



1/3 



proton 
3p 



3/2 



5/2 , ■ 



- - 2f. 



(114) 



li 



13/2 



2f™ 



_2Z2^ ■ 



^ \_ SLy5 +TF J \_ T31 +TF _ 



proton 
3p 



3/2 



5/2^ " 



(114). 



13/2 



lh„ 



IT. 



Ill 



9/2 - 



C 

•1— I 

I 

u 



2 - Ik 



-6 - 



neutron 



'17/2. ■ 



(184) 



2h. 



11/2 



4s 



3d 



3d. 



3/2 



8 -SLy5 +TF- 



neutron 
Ik. 



17/2^ ■ 



(184) 



2h . - 



11/2 



4s 



1/2 



3ff 



3d. 



■3/2 



T31 +TF 



FIG. 7: Single-particle spectra of ^^*114 near the Fermi sur- 
face. The calculations are performed with (SLy5-|-TF and 
T31-I-TF) and without (SLy5 and T31) the tensor interaction 
in the spin-orbit potential. 



V. SUMMARY 

A detailed investigation for the effect of the tensor force 
on the pseudospin energy splitting have been carried out 
taking Sn isotope chain as an example in the framework 
of the Skyrme-Hartree-Fock approach with the param- 
eter sets SLy5-HTF and T31-f-TF combined with BCS 
method, where the sign of ar and /3t pla-ys a crucial role. 
The contributions of the tensor force to the pseudospin 
symmetry are significant since all pseudospin doublets 
except (npi/2, np3/2) are composed of a j> orbit and a 
j< orbit which are affected in an opposite way by the ten- 
sor force. For {iypi/2, 117)3/2), only np3/2 orbit is shifted 
due to the tensor force, and hence the tensor force also 
affects the pseudospin splitting. Thus a lot of phenom- 
ena being related to the pseudospin symmetry should be 
accordingly affected. In addition, the effect of the tensor 
force on the shell correction energy is calculated taking a 
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magic nucleus ""^^^Sn and a superheavy nucleus 298114 as 
examples. It was shown that the influence of the tensor 
force on the shell correction energy is considerable. For 
2^^114, the absolute values of calculated proton shell cor- 
rection energy is enlarged due to the inclusion of the ten- 
sor component, which is consistent with the result that 
the shell gap at = 114 is obviously enhanced with the 
present parameter sets. The modification of the shell cor- 
rection energies by the tensor component are important 
for superheavy nuclei, which are related to their stabil- 
ity. All the conclusions here actually originate from the 
shifts of the single-particle levels (i.e. modified spin-orbit 
splittings) on account of the presence of the tensor force. 
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